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ABSTRACT 

The issue whether Moreton waves are flare-ignited or CME-driven, or a com- 
bination of both, is still a matter of debate. We develop an analytical model 
describing the evolution of a large-amplitude coronal wave emitted by the ex- 
pansion of a circular source surface in order to mimic the evolution of a Moreton 
wave. The model results are confronted with observations of a strong Moreton 
wave observed in association with the X3.8/3B flare/CME event from January 
17, 2005. Using different input parameters for the expansion of the source re- 
gion, either derived from the real CME observations (assuming that the upward 
moving CME drives the wave), or synthetically generated scenarios (expanding 
flare region, lateral expansion of the CME flanks), we calculate the kinematics of 
the associated Moreton wave signature. Those model input parameters are deter- 
mined which flt the observed Moreton wave kinematics best. Using the measured 
kinematics of the upward moving CME as the model input, we are not able to 
reproduce the observed Moreton wave kinematics. The observations of the More- 
ton wave can be reproduced only by applying a strong and impulsive acceleration 
for the source region expansion acting in a piston mechanism scenario. Based on 
these results we propose that the expansion of the flaring region or the lateral 
expansion of the CME flanks is more likely the driver of the Moreton wave than 
the upward moving CME front. 

Subject headings: shock waves — Sun: corona — Sun: flares 
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Introduction 



Solar flares and coronal mass ejections (CMEs) are explosive processes that are ab le 



to ge nerate large-scale wave-like disturbances in the solar atmosphere (e.g. 



Warmuth 



20071 ). Signatures of such disturban ces were first im aged in the hydrogen Ha spectral 



line and called Moreton waves after 



MoretonI ( 1960 



see also Moreton & Ramsey, 1960). 



Typically, Moreton waves appear as propagating dark and bright fronts in Ha filtergrams 
and dopplergrams, respectively, which can be attributed to a compression and relaxation 



of the chromosphe ric plasma. The disturbance propagates wit 



1000 km s"^ 



Veronig et al. 



(e^ 



Moreton &: Ramsey 



196G 



Zhang et al 



1 a speed in the order o f 



2001 



Warmuth et al 



2004a 



20061 ). which led to the conclusion that such a phenomenon cannot be of 



chromospheric origin, but is the surface track of a coro nal disturbanc e compressing the 



underlying chromosphere (sweeping- skirt hypothesis; see 



Uchida 



19681 ). M oreton waves are 



genera lly observed to be closely associated with the fiare impulsive phase (jWarmuth et al. 



2004a 


), whic 


1 often coincides a 


so wit 


1 the accelera 


tion p 


Zhang et al. 


2001; 


Vrsnak et al. 


2004; 


Maricic et al. 


2007; 



lase of the associate d CME (cf 



Temmer et al. 



20081). 



Moreton waves are observed to propagate perpendicular to the magnetic field, and 
the initial magnetosonic Mach numbers are estimated to lie in the range of M, 
suggesting that they are at 



2004 



Warmuth et al. 



range of M^s ~' 


.4-4, 


(Narukage et al. 


2002 



2004bl ). In their late propagation phase the wave perturbations 
undergo a broadening, weakening, and deceleration until Mms ~1 is reached. These results 
indicate that Moreton waves are a consequence of shocks formed from large amplitude-waves 



that decay to ordinary fast magnetosonic waves, which is in line with the fiare initiatec 



"blast wave" scena.rio (e . g., 



2002 



Vrsnak et al. 



2002 



iWarmuth et al 



Hudson et al. 



2003 



2001 



Khan fc Aurass 



Narukage et al. 



200^ 



Narukage et al. 



20041 ). Further evidence for 



the close association to shocks is the quasi-simultaneous appearance of Moreton waves 
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and radio type 


I bursts, which are one of t 


Khan & Aurass 


2002 




J'ohiolainen et aL 


2001 


2005; 


Vrsnak & Oliver 


2008 


)• 



2008 



Warmuth et aL 



2004b : 



Vrsnak et al 



Wave-hke disturbances were for the first time imaged directly in the corona by the EIT 



instrument aboard the So 



waves (IMoses et al 



1997 



ar and Heliospheric O bservatory (SoHO), thereafter called EIT- 



Thompson et al 



manifestation of the Moreton wave (IThompson et al 



19981) ■ They we re considered to be the coronal 



19991 ). but statistical studies revealed 



discrepancies in their velocities. E. 



than Moreton waves (IKlassen et al 



T waves were found to be two to three times slower 



2000l ). Today, their relation to Moreton waves and the 



generation mechanism of 



Wills-Davev &: Thompson 



IT waves is very much deba t ed (e.g.. 



1999 



Ohen &: Shibata 



2000 



Delannee fc Aulaniei 



Biesecker et al 



2002 



1999 



Oliver et al 



2004 



2007 



Veronig et al 



20081) 



Warmuth et al. 2004b: Oliver et al. 2005; Vrsnak et al. 2005; Ohen 2006; Attrill et al 



In the present paper, we solely focus on Moreton waves, which are generally accepted 
to be a chromospheric response to coronal shock waves. In particular, we study their 
generation mechanism and address the issue whether they are fiare-ignited or OME-driven, 
or a combination of both, which is still a matter of debate. To this aim, we developed a 
simple analytical model which describes the launch and propagation of Moreton waves. 
(Note that the presented model does not intend to evaluate generation mechanisms which 
may cause EIT waves.) We use for the model different input parameters acting as source 
that drives wave, first derived from OME observations (assuming that the upward moving 
OME drives the wave), and second using synthetically generated scenarios (to emulate 
alternative driving mechanisms). By confronting the results derived from the model with 
observations we aim to find constraints on the possible drivers of the wave. For this we use 
the outstanding observations of the Moreton wave associated with the X3.8/3B fiare-OME 
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event from January 17, 2005. We emphasize that the event was characterized by a very 
distinct a nd fast Moreton signatu re, indicating that it was caused by a coronal fast-mode 



shock (c.f. 



Warmuth et al. 



2004a) 



The observations of the Moreton wave and associated CME and flare under study are 
presented in Sect. 2. The model is described in Sect. [31 The results are given in Sect. 4. 
Discussions on the results, constrains of model input parameters by observations, and final 
conclusions are presented in Sect. 5. 



Observations 



Associated with the January 17, 2005 3B/X3.8 flare event, a fast Moreton wave starting 
at ~09:44 UT was observed with high time cadence (<1 min) in full-disk Ha filtergrams at 
Kanzelhohe Observatory. The wave propagated at a mean velocity of 930 km s~^ up to a 
distance of 500 Mm from its source location (f or more details on th e wave measurements 



and its propagation characteristics we refer to iVeronig et al. 



20061 ). The flare and its 



associated coronal mass ejection (CME) occurred at [N15,W25]. From this region actually 

bcusing on the 



Brueckner et al. 



two fast CMEs were launched within short time and in our study we are 
secon d event. The Large Angle and Spectrometric Coronagraph (LASCO; 
19951 ) instrument C2 aboard the Solar and Heliospheric Observatory (SoHO) imaged the 
first CME at 09:30 UT and the second CME at 09:54 UT. The linear plane-of-sky speed of 



the first CME was of ~2100 km s"^ and of the second CM 



with LASCO C2 and C3 (LASCO catalogue; lYashiro et al. 



of ~2500 km s ^ as observed 



20041 ). The study is performed 



over the time range 09:30-09:54 UT, hence, in the interval of interest we assume no impact 
on the CME kinematics due to the possible merging process with the previous event. 



The early CM: 



Hill et al. 



evolution could be observed with the G0ES12 Soft X-ray Imager (SXI; 



20051 ). Rising CME loops could be identified in 9 SXI frames with high time 
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cadence (~2-4 min; see lTemmer et al.ll2008l ) 



After co-aligning the GOES/SXI and Ha observations, the distances of the CME 
leading edge as well as the Moreton wave fronts were measured using as null-point the wave 



"radiant point' 
(for details see 



, which was derive d from circular fits to the earliest observed wavefronts 



Veronig et al. 



20061 ). From running ratio SXI images the height-time profile 



of the erupting CME structure was measured. 

In Fig. [1] we show the propagation of the Moreton wave together with the associated 
CME during its initial phase up to ~1 and the flare hard X-ray ( HXR) flux meas ured 
with the Ramaty High- Energy Solar Spectroscopic Imager (RHESSI; iLin et al.ll2002l ) in 
the non-thermal energy range 30-100 keV. From the second derivative of the height-time 
measurements we determined the onset of the CME fast acceleration phase, i.e. the 
launch time of the CME, at ~09:40-09:42 UT. The b ack-extrapolated Mo reton wave as 



well as the first HXR burst started at ~9:42 UT (see 



Veronig et al 



20061 ). The CME 



acceleration reached its peak of 4.4 ± 0.3 km s ^ at ~09:46 UT, and ends at ' ^10:06 UT (c f. 



jemmeretal 



(120071 ) and 



20081) . For the fu ll CME kinematics up to 30 R© we refer to 



Vrsnak et al. 



Temmer et al. 



tooi ). 



A composite dynamic radio spectrum for that day over the frequency range 600 MHz- 
20 kHz combining Artemis, DAM and WAVES measurements can be found under 



http : / /secchirh . obspm . f r/ select . php| The radio signatures show a rather complex 
situation most probably due to the launch of two CMEs for which a detailed study is 



given by 



Bouratzis et al. 



(120091 ). Associated with the event under study was a metric type 



II radio burst at 09:43-09:46 UT reported from San Vito, Italy (SVTO; spectral range 



"^Note that in 



Temmer et al. 



(120081 ) the Sun-center was used as null-point for the distance 



measurements of the CME. 
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70-25 MHz) and also from Learmonth, Australia at 09:44-09:47 UT (LEAR; spectral range 
65-25 MHz) as reported from the Solar Geophysical Data (SGD) under Solar Radio Spectral 



Observations ( ftp : //f tp . ngdc . noaa . gov/STP/SGD/ [ ) . Both stations report shock velocities 
of 1500 km s~^ using a one-fold Newkirk model which is consistent with an MHD shock 
moving through the solar corona. A group of type HI bursts occurred 09:41-09:47 UT, 
matching the main RHESSI peak. 

In Fig. [2] the distance-time and velocity-time profile of the observed Moreton wave is 
shown. The profile shows an increase in velocity with an initial speed of 400 km s^^ until it 
reaches a maximum speed of 1100 km s~^ at ~09:47 UT, afterwards the velocity decreases. 
This temporal behavior can be interpreted as nonlinear evolution of the wavefront. First, 
the wavefront steepens until a discontinuity appears, i.e. the shock formation starts. Then 
follows a phase of shock amplitude gro wth, which is refiecte d in shock acceleration and 



intensification (see Figures 4 and 5 in 



Vrsnak &: Lulic 



2000|). Finally, after the shock 



am plitude attains its maximum, the wave gradually decays to an ordinary fast-mode wave 
(cf. Izic et ahlboosl ). 



Fig. [3] shows the derived Moreton wave fronts with respect to the photospheric magnetic 
field. The first wave appearance is clearly located outside the active region. Since the wave 
propagated well outside the active region, Alfven speeds for the corona can be considered 



to lie in the range of 300-600 km s ^ (e.g. 



Narukage et al 



2002 



Warmuth fc Mann 



20051). 



The high velocity of the wave within a low Alfven speed environment as well as the 
associated metric type I I radio burst suggest that the wave is at least initially shocked (e.g. 



Gopalswamy et al.l 



1998 



Klassen et al. 



19991). 



The main criteria derived from the observations which our model results have to meet 
are 1) general kinematics of the wave, 2) velocity evolution and 3) timing of the shock 
formation. 
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3. The model 

We would like to emphasize that the following analytical model is kept as simple as 
possible and can thus only reproduce the general characteristics of the propagation of the 
disturbance. The model will simulate the Moreton wave by applying a driver which is a 
circular source region that may expand and move translatory at the same time. 



hree types of source expansion are applied following the terminology by IVrsnak 
( I2OO5I ): 1) The radius of the source is kept constant, i.e. there is no expansion of the source 
in time during its upward motion. Accordingly, plasma can flow behind the driver and 
the source acts as blunt-body driving a bow-shock. 2) The source radius expands with a 
constant radius-to-height ratio, r{t)/h{t), acting as a combined bow-shock/piston driver. 3) 
The source expands only in lateral direction without upward motion and plasma can not 
flow behind the contact surface, according to which the driver acts as piston mechanism. 

Our first intention is to investigate whether the Moreton wave could be produced by 
the upward moving CME, using the height-time measurements derived from the CME 
observations as input for the expanding source. We consider this model input for scenarios 
where the source acts as bow-shock and combined bow/piston driver for the wave (different 
strengths and proportions between the upward motion and lateral expansion of the driver 
are applied). Our second intention is to emulate an expanding flare region or the lateral 
expansion of the CME flanks for which we use synthetic expansion profiles. Such kind 
of model input is considered for a source that acts as piston driver mechanism for the 
wave. The results from the model will be compared to the kinematics of the January 17, 
2005 Moreton wave to estimate what kind of source expansion reproduces the general 
characteristics of the observed wave kinematics best. 

We suppose that the source accelerates to a high velocity, which causes a large 
amplitude coronal disturbance that is capable of compressing the underlying chromosphere 
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to produce the Moreton wave. The term large-amphtude waves should emphasize that 



the wave evolution can not be described through li nearized equat ions. For more detai 



VrsnakI (120051 ) and 



Warmuthl (120071 ). 



on the terminology of large scale waves we refer to 
In the case of a large amplitude wave, the rest frame velocity w of a given wavefront 
element (hereinafter called "signal") depends on two quantities. First, it depends on the 
local magnetosonic speed v^g, which is larger than in the unperturbed plasma due to the 
plasma compression, and is thus related to the perturbation amplitude. Second, it must 
be taken into account that a given signal propagates through a moving plasma, since the 
plasma flow velocity u associated with the perturbation amplitude is not negligible (see 
Fig. Bb). Consequently, the rest frame velocity of the signal equals to w = f ms + u (see 



Landau fc Lifshitz 



19871 ). i.e., elements of larger amplitude propagate faster. Due to the 



nonlinear evolution of the wave front, its profile steepens an d after a certain time/distance a 



discontinuity forms, marking the onset of shock formation (iLandau fc Lifshitziil987 



1995 



Vrsnak &: Lulic 



200G 



Zicetal 



20081). 



Mann 



Generally, the dependence of Vms on the perturbation amplitude cannot be expressed 
straightforwardly. However, in the case of a low plasma-to-magnetic pressure ratio (3 which 
is assumed here, the relationship simplifies, since the Alfven velocity va is much larger 
than the sound speed, and under the frozen-in condition in the case of perpendicular wave 
propagation, th e plasma density p is p roportional to the magnetic field strength B, i.e. 



Va oc y^p. 



Vrsnak fc Lulid (120001 ) have shown that in such a situation the relationship 



between the local propagation speed and the amplitude becomes very simple: the local 
value of Va can be expressed as va = vao + u/2, where f ao is the local Alfven velocity in 
the unperturbed plasma. Bearing in mind that w = va + u, one finally finds that the wave 



element propagates at the rest frame speed w = vao + Su/2 ( jVrsnak fc Lulic 



2000). 



Since the phase velocity of the signal depends on its amplitude u and the ambient 
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Alfven velocity fAo, the evolution of the wavefront depends on the spatial distribution of 
f AO and the evolution of the amplitude. The simplest possible situation is propagation of 
the wave in a medium where f ao is uniform. In such a case, the phase velocity changes only 
due to the amplitude evolution, which is governed by the energy conservation. For example, 
in the case of a spherically symmetric source, creating a spherically symmetric wavefront 
(Fig. IDd), the amplitude is inversely proportional to the di stance d, i.e., decreases as d"^, 



whereas in the cylindrical symmetry it decreases as d ( Landau fc Lifshitz 



19871) • Note 



that in the case of freely-propagating shock waves (blasts), the amplitude decreases also 
because the leading edge of the perturbation (having the highest velocity) propagates faster 
than the low-amplitude segments in the trailing edge. This causes perturbation profile 



broadening, which must be compensated by an amplitude decrease (iLandau fc Lifshitz 
1987)@ 



Of course, in the solar corona the Alfven velocity is far from being uniform. Even if 
the coronal structural inhomogeneities ar e neglected, it chaiiges w ith height and depends 



on the distance from active regions (e.g., IWarmuth fc MannI 120051 ). In such a situation 



where the spatial distribution of fAo is generally unknown, one has to investigate the 
wavefront kinematics by calculating the amplitude evolution for various reasonable spatial 
distributions of vao- However, instead of this, we apply an analogous procedure, where we 
take f AO uniform, and describe the signal amplitude and the phase-velocity evolution by 
different functional forms. In other words, instead of presuming a function that describes 
the change of vao with distance d from the wave source, we directly presume a function that 



^Note that in a medium where the Alfven velocity decreases steeply enough with the 
distance, the leading edge might be slower than the trailing edge. In such a case, the 
wavefront slows down, whereas the amplitude increases. 
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describes the wave evolution. In particular, we use the power-law function 

f{d) = d-^ (1) 

and exponential function 

f{d)=e-'/^. (2) 

Applying different decay lengths (denoted in the power-law function by a and in the 
exponential by p) we can reproduce a weak or strong attenuation of the signal. Note that 
/ = 1 would represent a plane wave without decay as achieved for p — > oo and a — > 0. On 
the other hand, large a or small p represents a strong attenuation. 

Beside the power-law and exponential function, we also employ as a kind of reference, 
the functions: 

f{d) = (3) 



and 



f{d) = ^ , (4) 



which describe the amplitude decrease of cylindrically and spherically symmetric sound 
waves, respectively. 

The initial amplitude of a given signal is determined by the velocity of the source 
surface v^. At the starting time to when the signal is launched, u{to) = Vs(to), since the 
flow velocity has to be equal to the contact-surface velocity. The geometry of the source 
is considered as a radially expanding surface of cylindrical (arcade expansion) or spherical 
shape (volume expansion) with a radius r{t) centered at the height h{t). Applying the 
Huygens-Fresnel principle, one finds that due to the presumed symmetry of the source and 
the presumed homogeneity of the ambient plasma, the wavefront elements are concentric 
with the source surface (cf. Figs. Sb, O and [S]). 



We follow the signals which are emitted continuously from the source surface for the 
time span until a certain time at each small time step At = tj — The distance x 
traveled by the signal from the time to when it was emitted, until the time tj, is calculated 
iteratively. Using the expression 



we obtain the distance from the source region center, d{ti) = r{to) + x{ti), where r(to) is 
the radius of the source surface at the time to; when the signal was emitted (Fig. |Dd). Note 
that x(to) = and c?(to) = "'"(^o); and Eq. Ohas to be integrated from to to tc- 

Considering the mimicked Moreton wave as the extension of the outermost signal 
measured at the solar surface (cf. arrows in Figs. [S]and[H]), we derive for each time step At 
the propagation of the wave as distance du{t)- Hereinafter, this outermost signal that is 
considered to mimic the Moreton wave, will be denoted as the ground track signal (GTS). 



In the following, distance-time plots and velocity profiles are shown for the propagated 
GTS resulting from our model. The results are confronted with the observed Moreton wave 
kinematics. Due to the huge spectrum of possibilities obtained by varying and combining 
the different model parameters, we will show here only representative model results, i.e. 
those which match the observational criteria of the Moreton wave best. The successful 
model will reproduce the general characteristics of the observed Moreton wave in terms of 1) 
kinematics, 2) velocity evolution (increasing velocity until ~09:47 UT followed by decreasing 
velocity), and 3) shock formation around the onset of the type II burst (~09:43 UT), i.e., 
before or close in time to the first appearance of the Moreton wave (~09:44 UT). 




4. Implementation and interpretation of the model 



- 13 - 



The wave-like disturbance that generates the Moreton wave is assumed to propagate 
approximately near the cor onal base. Under this ass umption, the value for vaq lies in the 



range of ~300-600 km s ^ (jWarmuth &: Mann 



20051 ). To ease the comparison between the 



model results and the observations (bearing in mind also other aspects of the CME/flare 
event) we use for the model the absolute time in UT. The parameter to varies around 
~09:42 UT which is close to the onset of the fast acceleration stage of the CME and the 
flare onset in Ha and HXRs. The paramet er tr is the tirae at which the Moreton wave 



was observed the last time (~9:54 UT; see 



Veronig et al 



20061 ). The time range t^-tc is 



subdivided into time steps At=10 s, i.e. each 10 seconds the position of the wavefront and 
the GTS is calculated. 



4.1. Model results based on observed CME kinematics 

In Fig. we give a snapshot of the propagated signals (circles) that were emitted 
during the upward motion (along the y-axis) of an expanding source. The kinematics 
for the upward moving source is taken from the CME observations, and the type of 
source expansion acts as a combined bow shock/piston driver for the emitted signals with 
r{t)/h{t) =0.2, i.e. source size is proportional to height at each time t. The decay of the 
signal is based on a cylindrical geometry of the source (see Equ. [3]). The first signals are 
emitted at to=9:41:52 UT when the CME had a height of /i(to)=105 Mm and an initial size 
of r(to)=21 Mm. The surrounding Alfven speed of the unperturbed plasma is chosen as 
fAo=400 km s~^. From to on, we follow the signals every 10 s, until they have reached a 
certain extension at tc=9:53:52 UT (Fig. [5]). Note that signals which are launched right after 
to have the longest time to evolve, signals launched close to tc the shortest. At 9:53:52 UT 
the CME has a height of h{tc)=1570 Mm and a size of r(tc)=314 Mm. The arrow in Fig. [5] 
indicates the propagated distance (iM(^c)=881 Mm of the GTS, i.e. the mimicked Moreton 



-14- 



wave at 9:53:52 UT. 

Fig. O shows the calculated GTS distance versus time using the observed CME 
kinematics as input for the upward moving source for two different types of the source 
expansion. The top panel of Fig. [6] is supposed to mimic a combination of a bow shock and 
piston driven scenario; the source was expanding during its upward motion self-similarly 
with a constant ratio of r{t)/h{t) = 0.6. The bottom panel of Fig. [6] supposes the source to 
act as a rigid-body driver, i.e. the radius was kept constant during its upward movement 
with r(t)=140 Mm, imitating a bow-shock scenario. 

The derived kinematics of the GTS show a distinct feature of a "knee" as indicated in 
the top panel of Fig. [61 The feature occurs when a later e mitted GTS passes th e preceding 



onqfl, i.e. the knee marks the time of the shock formation ( Vrsnak fc Lulic 



2000). 



From Fig. [6] it can be seen that the first phase of the observed Moreton wave could 
be partly mimicked but not its later evolution. The knee, which represents the time of the 
shock formation, occurs ~4-6 minutes after the first Moreton wave front was observed. 
In Fig. [7] the according velocity profiles are plotted for the scenarios presented in Fig. [6l 
For both scenarios, CME acting as combined bow/piston and bow driver, the GTS is 
of decreasing velocity until ~09:47 UT and the velocity of the GTS at ~09:51 UT (last 
observational data point) is about 1.5 times as high as for the observed Moreton wave. 
Hence, the CME is a too fast driver which generates a too fast GTS at large distances. 
Although various kinds of parameter values were applied, it was not possible to reproduce 
the general observational characteristics of the Moreton wave. From this we conclude that, 
using a fast upward moving driver for the model, like the observed CME, all generated 



■^In the specific case of our model the overtaking GTS was launched when the source 
speed changed from subsonic to supersonic. 
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GTS profiles reveal 1) increasing velocity after ~09:47 UT and 2) a shock formation 
several minutes after the first observed front of the Moreton wave (cf. Fig. [2]), which is not 
consistent with the observations. 

4.2. Model results based on a synthetic kinematical profile of the source 

From the calculated GTS kinematics using real CME observations, it became clear that 
the radially upward moving CME, imitating a bow or combined bow/piston scenario, cannot 
reproduce the observed Moreton wave characteristics. In order to investigate alternative 
driving mechanisms, we use as input parameters a synthetic kinematics of an expanding 
source acting as piston mechanism. As simplest approach, we assume that during the radial 
expansion the center of the source is fixed at the surface, i.e. h{t)=0, in order to imitate 
a spherical or cylindrical piston. The synthetic kinematics consists of an acceleration 
phase ta of constant acceleration a, until a certain velocity is reached by which the source 
expands further. This enables us to study the signal evolution emitted from very differently 
expanding driving sources, ranging from sudden impulsively to gradually accelerating. 

In Figs. M and [9] a relatively gradual expansion of a spherical piston is represented. 
We use as input an initial source size of r (to) =140 Mm accelerating over a time span of 
ta=AOO s with a=2.8 km s^^ (final velocity 1120 km s^^). The arrow in Fig. [8] indicates 
the propagated distance (iM(^c) of the GTS, i.e. the mimicked Moreton wave. The shock 
formation time was obtained at 09:48:32 UT, hence, several minutes after the first Moreton 
wave front was observed. The corresponding velocity profile (shown in Fig. [11] as dashed 
line) reveals an increase of velocity of the GTS in the late propagation phase after 09:47 UT, 
although a strong decay (exponential) was applied to the GTS its kinematics. Similar to 
what we obtained applying the observed CME kinematics such an acceleration behavior of 
the source cannot mimic the observed Moreton wave. To compensate for the delayed timing 
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of the shock formation, a shorter and more impulsive acceleration of the source expansion 
would be required to reproduce adequately the Moreton wave propagation. 

The top panel of Fig. [TO] shows the expansion of a smaller source of r(to)=110 Mm of a 
shorter and stronger acceleration (ta=160 s; a=4.8 km s~^) in comparison to the previous 
scenario. The calculated GTS from this case shows a very good match with the observed 
Moreton wave kinematics as well as its velocity profile (dotted line in Fig. [TTl) . The timing 
of the shock formation at 9:44:32 UT is close to the first detected Moreton wave front 
(~9:44:30 UT). Since after the shock formation the GTS propagates faster than the later 
emitted signals we assume that the source is acting only temporarily as piston. The time 
range during which the wavefront evolves independently from the driver is indicated as 
dashed gray line in Fig. [101 From this we derive, the source surface would need to expand 
from the initial size of 110 Mm up to 170 Mm to mimic the resulting Moreton wave (solid 
gray line in Fig. [TOj) . The initial source size of ~110 Mm would roughly correspond to the 
diameter of the active region (cf. Fig. [3]). A further scenario is presented in the bottom 
panel of Fig. [TO] with source parameters comprising an initial size of r(to)=50 Mm and a 
very impulsive expansion (short and strong acceleration) of ta=80 s and a=8 km s^^. The 
synthetic kinematics of the calculated GTS matches the observed Moreton wave reasonably 
and the shock formation for this scenario takes place at 9:42:52 UT. The source surface, 
acting as a temporary piston, would need to expand from its initial size of 50 Mm up to 
75 Mm. Considering the velocity profile (dashed-dotted line in Fig. [TTl) the GTS reaches its 
peak velocity before 09:47 UT, however, decreases very rapidly. Compared to the earlier 
scenario (source parameters: r(to) = 110 Mm; ta=lQO s; a=4.8 km s~^; marked with the 
dotted line in Fig. [TTl) ^^e match is worse, however, still reasonable within the limits of such 
a simple model. 

In Fig. [TH we show the velocity profiles from the simulated wave kinematics as given in 
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Figs. [9] and [T0|, and compare them to the velocity profile derived from the observed Moreton 
wave (solid line). We obtain the best match for a wave which is assumed to be driven by 
a shortly and strongly accelerating source (dotted line); a more impulsive expansion of 
the source would generate a profile of comparable velocity at the last point of observation 
close to 09:51 UT, but peaks earlier (dashed-dotted line). Such source behavior could be 
interpreted as the expanding flanks of a CME or the volume expansion of a flare. On the 
other hand, a weak and long acceleration similar to the upward moving CME (dashed line) 
reveals substantial inconsistencies to the observed wave profile (late peak, final velocity too 
high). 



Discussion and Conclusion 



The analytical model presented here is based on tracing the evolution of a large 
amplitude wave. This is justified since Moreton waves are caused by a strong compression 
of the chromosphere (otherwise the wave would not be seen in Ha). There are several 
unknown factors whose implementation would be beyond the scope of this model. For 
example, we considered a homogeneous corona where the density and the Alfven velocity do 
not change, neither in the vertical nor in the horizontal direction, taking vao in the range 
300-600 km s~^. Recent observational studies showed that the magnetosonic speed v^as (we 
assume v^s ~ "^^ao) can drop down to a local minimum of 300-500 km s^^ around the height 



2 Rq but t hen rises steadily up to a local maximum o 



3 and 4 R© (IMann et al. 



2003 



Warmuth fc Mann 



20051 ). 



1000 km s ^ a t a he ight between 



Vrsnak et al. 



(12004h obtained 



from observations of type II bursts that on average the magnetosonic speed attains a local 
minimum of f ms ~400 km s~^ around 3 R© and a broad local maximum of v^s ~500 km s^^ 
in the range of 4-6 Rq. Besides, the previous CME event which started about 40 min 



earlier (LASCO catalogue; 



Yashiro et al. 



20041 ) from the same active region might affect the 
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actual value of the Alfven velocity too. 



Furthermore, we did not take into account the accurate relation between the plasma 
flow and source velocity u, i.e. the CME velocity, but simply used a one-to-one relation. We 
approximated u by the CME speed which is appropriate concerning the upper part of the 
moving and expanding CME but does not hold for the lateral direction, i.e., from which the 
GTS kinematics is determined. We tried to account for this by reducing the CME speed 
by ~60%, thus maintaining the CME kinematical profile as model input but with a lower 
speed. However, also that option did not result in a better match between the generated 
GTS and the observed Moreton wave. 

An important factor for the derived model results is the decay factor used to attenuate 
the signal. Since in the corona the distribution of density p(r), magnetic field B{r), and 
Alfven speed VAo{r) are unknown, we use different "decay functions" (see Equ. [T]-|l]). It 
had to satisfy two criteria: it should be strong enough to decelerate the signal in its 
late propagation phase but should not, due to its strength, delay the timing of the shock 
formation. We used geometry dependent factors adapted from sound wa ves (cylindrical 



Zic et al. 



20081 ). 



and spherical), i.e., without implementing a magnetic field (for details see 
Formal decay factors, like power-law and exponential functions, were used to put the decay 
to the limits either having no attenuation o r very strong attenuation and to account for the 



unknown distribution of Vy, 



Pagano et al 



(120071 ) investigated the role of magnetic fields 



for an expanding and upward traveling CME and showed that a spherical cloud without 
a magnetic field drives a wave that propaga tes to longer distances than that with a weak 



open field (see Fig. 7 in 



Pagano et al. 



20071 ). This implies that the presence of a magnetic 
field would result in a stronger signal decay than obtained from our simple approaches. 
Since we were not able to reproduce the wave using the limits for the decay factor (strong 
versus no attenuation), we suppose that even utilizing more sophisticated decay factors, the 



- 19 - 



disturbance generated by the CME forehead would not be able to reproduce the observed 
Moreton wave. 

Using the observed CME kinematics as input parameters the model could not 
reproduce the general characteristics of the observed Moreton wave. The timing of the 
shock formation ("knee") was not appropriate but occurred later than the first observed 
Moreton wave front. The velocity profile was not conform and the final velocity was too high 
in comparison to the observed Moreton wave. By varying the initial source size as well as 
the behavior during the source evolution (bow, piston or combined bow/piston driver), the 
GTS kinematics was shifted to a larger or sm aller propagatio n distance, however, the shock 



formation always appeared too late (see also 



ZicetaL 



20081 ). Similar results are obtained 



by applying different start times for the signal to and different local Alfven velocities wao- 
Thus, experimenting with all these different parameters demonstrated that the Moreton 
wave could not be reproduced when taking the kinematics of the radial outward movement 
of the CME as input for the model. 

This finally pushed us to use synthetic kinematics in order to imitate other possible 
drivers for the signal. So far it was clearly derived from the model that the source expansion 
needs to be more impulsive (early shock formation). For synthetic kinematics of stronger 
and shorter acceleration of the source surface expansion (3-D piston type) we found a good 
match between the model generated signal and the observed Moreton wave. The timing 
of the shock formation is, when using these kinematical profiles, in good agreement with 
the appearance of the first Moreton wave front. Using an exponential attenuation factor 
(see Equ. [2]) with short signal decay lengths the best match to the observed Moreton wave 
could be found. On average the Alfven Mach number Ma from such synthetic kinematics 
are within the range of Ma ~ 1.5-3 which agrees with observed A lfven Mach numbers for 



Moreton waves (e.g. 



Narukage et al 



2002 



Warmuth et al. 



2004al ). The initial source size 
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and its expansion dimension that is necessary to mimic the observed Moreton wave can be 
interpreted as the laterally expanding CME flanks or the volume expansion of the flare. 



Pomoell et al. 



(120081 ) concluded from a 2D magnetohydrodynamic simulations that 
for the driver of a Moreton wave a high acceleration during a short time interval is 
necessary. This was interpreted to require a strong lateral expansion, either lif t-off of an 



Zicetal 



over-p ressured flux rope or thermal explosion-kind of energy release. Likewise, 
( 120081 ) obtained from a 3D analytical model that a short acceleration phase up to high 
velocities (~1000 km s^^) within a low Alfven velocity environment is necessary to create 
a shock that is capable of causing type II bursts in the dm/m wavelength range and Ha 
Moreton waves. 

Concluding, for the January 17, 2005 event under study it is unlikely that the bow 
shock of the CME generated the observed Moreton wave. The CME is a too gradually 
accelerating source in the lift-off phase and a too fast one in the later evolution phase to 
cause the observed Moreton wave kinematics. An impulsively accelerated expansion of a 
source surface acting as a temporary piston would be a more appropriate mechanism to 
generate the observed Moreton wave. Possible driving mechanisms would be the laterally 
expanding CME flanks or the impulsive volume expansion of the flare. The latter scenario 
would be in accordance with the flare initiated "blast wave" scenario proposed fro m 



observationa 



Vrsnak et al 



resu 



Warmuth et al 



Chen et al 



2001 



2004a 



(l2002h who 



ts for the kinematics of Moreton waves (see 
20021 ) but in contrast to the numerical model by 
claimed that Moreton waves correspond to the piston-driven shock over the CME. For 
the future it would be important to have more such complete data sets including both, 
observations from the early CME evolution (upward moving front as well as expanding 
flanks) and detailed observations of Moreton waves, in order to validate the presented 
results. 



- 21 - 



M.T. is supported by the Austrian Fonds zur Forderung der wissenschaftlichen 
Forschung (FWF Erwin-Schrodinger grant J2512-N02). T.Z. and B.V. acknowledge 
funding by the Croatian Ministry of Science, Education and Sports under the project 
007-0000000-1362. A.V. gratefully acknowledges the Austrian Fonds zur Forderung der 
wissenschaftlichen Forschung (P20867-N16). 



-22- 
REFERENCES 

Attrill, G. D. R., Harra, L. K., van Driel-Gesztelyi, L., & Demoulin, P. 2007, ApJ, 656, 
LlOl 

Biesecker, D. A., Myers, D. C, Thompson, B. J., Hammer, D. M., & Vourlidas, A. 2002, 
ApJ, 569, 1009 

Bouratzis, C, Preka-Papadema, P., Moussas, X., Alissandrakis, C, & Hillaris, A. 2009, 

Advances in Space Research, 43, 605 

Brueckner, G. E., et aL 1995, Sol Phys., 162, 357 

Chen, P. F. 2006, ApJ, 641, L153 

Chen, P. F., & Shibata, K. 2000, ApJ, 545, 524 

Chen, P. F., Wu, S. T., Shibata, K., & Fang, C. 2002, ApJ, 572, L99 

Chver, E. W., Laurenza, M., Storini, M., & Thompson, B. J. 2005, ApJ, 631, 604 

Chver, E. W., Nitta, N. V., Thompson, B. J., & Zhang, J. 2004, SoL Phys., 225, 105 

Delannee, C, & Aulanier, G. 1999, SoL Phys., 190, 107 

Gopalswamy, N., et aL 1998, J. Geophys. Res., 103, 307 

Hill, S. M., Pizzo, V. J., Balch, C. C, Biesecker, D. A., & et al. 2005, Sol. Phys., 226, 255 

Hudson, H. S., Khan, J. I., Lemen, J. R., Nitta, N. V., & Uchida, Y. 2003, Sol. Phys., 212, 
121 

Khan, J. I., & Aurass, H. 2002, A&A, 383, 1018 

Klassen, A., Aurass, H., Klein, K.-L., Hofmann, A., & Mann, G. 1999, A&A, 343, 287 



-23- 

Klassen, A., Aurass, H., Mann, G., & Thompson, B. J. 2000, A&AS, 141, 357 
Landau, L. D., & Lifshitz, E. M. 1987, Fluid Mechanics (2nd edition. New York) 
Lin, R. P., Dennis, B. R., Hurford, G. J., & et aL 2002, SoL Phys., 210, 3 
Mann, G. 1995, J. Plasma Phys., 53, 109 

Mann, G., Klassen, A., Aurass, H., & Classen, H.-T. 2003, A&A, 400, 329 

Maricic, D., Vrsnak, B., Stanger, A. L., Veronig, A. M., Temmer, M., & Rosa, D. 2007, 
Sol. Phys., 241, 99 

Moreton, G. E. 1960, AJ, 65, 494 

Moreton, G. E., & Ramsey, H. E. 1960, PASP, 72, 357 

Moses, D., et al. 1997, Sol. Phys., 175, 571 

Narukage, N., Eto, S., Kadota, M., Kitai, R., Kurokawa, H., & Shibata, K. 2004, in lAU 
Symposium, Vol. 223, Multi- Wavelength Investigations of Solar Activity, ed. A. V. 
Stepanov, E. E. Benevolenskaya, & A. G. Kosovichev, 367-370 

Narukage, N., Hudson, H. S., Morimoto, T., Akiyama, S., Kitai, R., Kurokawa, H., & 
Shibata, K. 2002, ApJ, 572, L109 

Pagano, P., Reale, P., Orlando, S., & Peres, G. 2007, A&A, 464, 753 

Pohjolainen, S., Hori, K., & Sakurai, T. 2008, Sol. Phys., 253, 291 

Pohjolainen, S., et al. 2001, ApJ, 556, 421 

Pomoell, J., Vainio, R., & Kissmann, R. 2008, Sol. Phys., 253, 249 



Temmer, M., Veronig, A. M., Vrsnak, B., Rybak, J., Gomory, P., Stoiser, S., & Maricic, D. 
2008, ApJ, 673, L95 

Thompson, B. J., et al. 1999, ApJ, 517, L151 

Thompson, B. J., Plunkett, S. P., Gurman, J. B., Newmark, J. S., St. Cyr, O. C, & 
Michels, D. J. 1998, Geophys. Res. Lett., 25, 2465 

Uchida, Y. 1968, Sol. Phys., 4, 30 

Veronig, A. M., Temmer, M., & Vrsnak, B. 2008, ApJ, 681, L113 

Veronig, A. M., Temmer, M., Vrsnak, B., & Thalmann, J. K. 2006, ApJ, 647, 1466 

Vrsnak, B. 2005, EOS Transactions, 86, 112 

Vrsnak, B., & CUver, E. W. 2008, Sol. Phys., 253, 215 

Vrsnak, B., & Lulic, S. 2000, Sol. Phys., 196, 157 

Vrsnak, B., Magdalenic, J., Temmer, M., Veronig, A., Warmuth, A., Mann, G., Aurass, H., 
& Otruba, W. 2005, ApJ, 625, L67 

Vrsnak, B., Magdalenic, J., & Zlobec, P. 2004, A&A, 413, 753 

Vrsnak, B., Maricic, D., Stanger, A. L., Veronig, A. M., Temmer, M., & Rosa, D. 2007, 
Sol. Phys., 241, 85 

Vrsnak, B., Warmuth, A., Brajsa, R., & Hanslmeier, A. 2002, A&A, 394, 299 

Warmuth, A. 2007, in Lecture Notes in Physics, Berhn Springer Verlag, Vol. 725, Lecture 
Notes in Physics, Berlin Springer Verlag, ed. K.-L. Klein & A. L. MacKinnon, 107 — h 

Warmuth, A., & Mann, G. 2005, A&A, 435, 1123 



-25- 

Warmuth, A., Vrsnak, B., Aurass, H., & Hanslmeier, A. 2001, ApJ, 560, L105 

Warmuth, A., Vrsnak, B., Magdalenic, J., Hanslmeier, A., & Otruba, W. 2004a, A&A, 418, 
1101 

Warmuth, A., Vrsnak, B., Magdalenic, J., Hanslmeier, A., & Otruba, W. 2004b, A&A, 418, 
1117 

Wills-Davey, M. J., & Thompson, B. J. 1999, Sol. Phys., 190, 467 

Yashiro, S., Gopalswamy, N., Michalek, G., St. Cyr, O. C, Plunkett, S. P., Rich, N. B., & 
Howard, R. A. 2004, Journal of Geophysical Research (Space Physics), 109, 7105 

Zhang, J., Dere, K. P., Howard, R. A., Kundu, M. R., & White, S. M. 2001, ApJ, 559, 452 

Zic, T., Vrsnak, B., Temmer, M., & Jacobs, C. 2008, Sol. Phys., 253, 237 



This manuscript was prepared with the AAS macros v5.2. 



- 26 - 



700 



600 - 



500 



400 - 



« 300 ^ 



200 



100 - 



09:32 



A SXI (CME sign.) 

+ Moreton wave 
RHESSI 30-100 keV 



Distance meosurements 
for CME and Moreton: 
center=[430, 295" 




09:36 09:40 09:44 09:48 

Start Time (17-Jan-05 09:32:00) 



09:52 



Fig. 1. — Evolution of the flaring process, CME take-off, and Moreton wave. Triangles 
indicate the measured CME height as derived from SXI images. Gray lines give the RHESSI 
flare hard X-ray flux for the energy band 30-100 keV. Plus signs give the kinematics of 
the observed Moreton wave front together with error bars. Around ~9:42 UT we observe 
the first Ha flare brightenings, the first HXR peak, as well as the CME and Moreton wave 
initiation (derived from back extrapolation). 
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Fig. 2. — Evolution of the Moreton wave kinematics (left scale) and velocity evolution (right 
scale). Plus signs with error bars show the Moreton wave front together with a 4*^ order 
polynomial fit (dashed line). The solid line shows the derivative of the polynomial fit, i.e. the 
velocity profile. The velocity curve can be divided into an early phase of increasing velocity 
(feariy) foUowcd by a later phase of decreasing velocity (fiate) with the inflection point at 
~09:47 UT (dashed vertical line). The gray horizontal bar indicates the local Alfven velocity 
fAo for the corona outside of active regions that lies in the range 300-600 km/s. 
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Fig. 3. — SOHO/MDI magnetogram scaled to a magnetic field strength of ±700 G. Solid 
white lines show the Moreton wave fronts observed in Ha images. The first wave front at 
09:44:30 UT is clearly located outside the area of strong magnetic fields. Plus sign indicates 
the wave radiant point derive d from a circular fit to the earliest observed wave front (for 



details see 



Veronig et al. 



20061). 
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Fig. 4. — a) Definition of tlie term "signal" : a given element (gray box) of the disturbance 
profile is referred to as signal. It is characterized by the propagation speed w and the 
associated plasma flow velocity u. b) The signal radius d at time tc, d{tc), is defined as the 
sum of the source size r at the time to when the signal was launched, and the propagated 
distance x of the signal up to the time tc (see Equ. [5]). 
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Fig. 5. — Black circles show the upward moving CME (observed CME kinematics) which 
is expanding during its upward movement from until with r{t)/h{t)=0.2 (combined 
bow/piston scenario). The cross indicates the height of the CME center at t^. Gray circles 
show the evolution of calculated signal radii launched from the circular source surface (CME) 
followed in steps of At=10 s. For time tc=9:53:52 UT the outermost signal measured at the 
surface h = 0, i.e. the mimicked chromospheric Moreton wave front, is marked with an 
arrow. 
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Fig. 6. — Calculated distance-time profile of the generated Moreton wave signal (solid black 
line) using observed CME kinematics as input for the model. The different parameter values 
used are specified in the legend. The solid gray line is the kinematics of the driver of the 
wave signal, i.e. in the present case the CME, the dashed gray line represents the stage when 
the driver is no more directly related to the evolution of the wavefront profile. The triangles 
indicate the CME front measurements from GOES/SXI observations. The kinematics of 
the observed Moreton wave front is plotted by plus signs with error bars. Top panel: a 
combined bow/piston driven scenario is applied using a ratio r{t)/h{t)=0.6 for the increase 
of the source size during its upward movement. Bottom panel: a bow shock driven scenario 
is applied, i.e. keeping the source size constant (r(t)=140 Mm) during its upward movement. 
The arrow pointing to the a;- axis indicates the launch time of the first signal tn. 



- 32 - 



2000 



£ 1 500 



o 
5 



1000 



500 



observed 
Moreton wove 

upword CME 
bow/ piston 
bow 




\ , , , \ , , 

09:44 09:48 

Start Time (17-jQn-05 09:43:00) 



09:52 



Fig. 7. — Velocity profiles derived from the generated wave kinematics shown in Fig. [6] 
using the observed kinematics of the upward moving CME as input. A combined bow/piston 
scenario (dashed-dotted line) and a bow shock scenario (dotted line) is applied for the driver. 
For comparison, the velocity profile of the observed Moreton wave is plotted as solid line. 
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Fig. 8. — Disturbance signals emitted from an expanding piston source using ta=400 s and 
a=2.8 km s~^. The arrow indicates the intersection of the outermost signal with the x-axis, 
i.e. the propagated way of the mimicked Moreton wave at 9:53:52 UT. The initial source size 
is r(to)=140 Mm. 
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Fig. 9. — Same as Fig. [6] but for a GTS derived from a synthetic kinematical profile of the 
source representing a spherical piston. The synthetic kinematical profile would be appropri- 
ate for a fast upward moving CME event. Different parameter values are specified in the 
legend. 



- 35 - 



L_ 

09:40 



v„: 400 km/s 
r^: 110 Mm 
t„: 160 s 
a: 4.8 km/s^ 
09:41:52 UT 



signal 

source 

otten: exponential 
decoy length: 550 Mm 




09:42 



09:44 09:46 09:48 

Start Time (17-Jan-05 09:40:00) 



09:50 




09:40 



400 km/s 
50 Mm 
80 s 

8.0 km/s^ 
09:41:32 UT 

signol 

source 

otten: exponentiol 
decay length; 600 Mm 



09:42 



09:44 09:46 09:48 

Start Time (17-Jan-05 09:40:00) 



09:50 



09:52 



Fig. 10. — Same as Fig. [9] but for synthetic kinematical profiles of different source sizes, 
acceleration times and strengths. Different parameter values are specified in the legends. 
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Fig. 11. — Velocity profiles derived from the simulated waves as shown in Figs. 191 and fTOl 
The dashed line represents the wave driven by a source region expansion of a=2.8 km s~^ and 
ta=400 s. The dashed-dotted line represents the source region expansion of a=8.0 km s"^, 
ta=80 s, whereas the dotted line shows the case a=4.8 km s~^ and ta=160 s. For comparison 
the velocity profile of the observed Moreton wave is plotted as solid line. 



